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Pseudo-rapidity distributions of high multiplicity JACEE events ( Ca-C, Si-AgBr 
and G27) and Texas Lone Star (TLS) event are analysed by the wavelet transform 
method. Using the Daubechies' waveltes, D4, Dq, Dg, and Dio, wavelet spectra of 
those events are calculated. Two JACEE events are compared with the simulation 
calculations. We discuss how to distinguish prominent peak from apparent peaks. 
The wavelet spectrum of Ca-C event seems to resemble that simulated with the 
Poisson random number distributions. The wavelet spectrum of TLS event shows 
, an oscillatory behavior. It is found that the wavelet spectra of JACEE G27 event 

\l ' are different from those of Arizona group. 
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^^ ■ 1 Introduction 

^ij . In high energy nucleus- nucleus (AA) collisions number density of secondary 

1-^ ' particles in the rapidity space becomes very high and studies of number den- 

r^. sity fluctuations in the rapidity spaceJa^ ^xpected to reveal new features of 

O ,' multiparticle production mechanisms EfaU'Ll The pseudo-rapidity distributions 

(D ! of the threcL^JACEE eventsi, (Ca-C, Si-AgBr and G27 ) and the Texaa^Lone 

i-^ ' Star (TLS)EI event are analysed by the wavelet transform methodLO'ElEEj. 

Any function (data) can be expanded into self-similar wavepackets in this 
scheme. Wavelet spectra of those four events are calculated and some of them 
. ■ . are compared with simulation calculationsu In § 2 the wavelet transform con- 

C^ ' cept is introduced and in § 3 the Daubechies' mother waveltes, D4, Dq, Ds, 

and -Dio, are prepared to be used in § 4. The wavelet spectra of the above 
mentioned four events are calculated and compared with the simulation calcu- 
lations in § 4. The final section is devoted to the concluding remarks. 

2 Wavelet transform 

Wavelets are constructed from dilation and translation of a scaling function 
4>{x), which is constructed by an iteration equation, 

N-l 

M^) = ^ Ck(t>.i-i{2x - k) i = l,2,---, (1) 

fc=0 
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from a primary scaling function (J)q (x) ( (/)o (x) = 1 for < 2; < 1 and (f>o (x) ~ 
otherwise). In Eq.dfl), N is a even number and Ck (/c = 0, 1, • • • , iV — 1) are 
constants. The iteration is continued until (jjii^) becomes indistinguishable 
from (f>i^i{x) and defines then 4'{x). 

The mother wavelet W{x) is given by 



Af-l 



^(^) = J2 (-l)''+'cAr_i_fc (b{2x + k~N+l), 



k=0 



and defines the j-th level wavelet (j = 0, 1, • • •) 

^j,kix) = 2iw{2^x-k), /c^O, l,•■•,2^'-l■ 



(2) 



(3) 



Coefficients Cfe (A; = 0, 1, • • • , iV — I) should be determined in such a way that 
wavelets and scaling function satisfy the following conditions, 

(f){x) 4>{x)dx = 1, / (p{x)'il}j_k{x)dx — 0, 

tpj,k{x)AAx)dx = SjrSks, I x''W{x)dx^0,k^0,---,N/2~l{4:) 

If N=4, we have 

l + ^/3 3 + V3 3-V3 1-^3 

Then arbitrary function f{x) defined in the region < x < 1 can be expanded 
as 

00 2^-1 

/ = aQ(j){x) + X! X! 0!j,k1p3,k{x)- 
j=0 fc=o 

We shall consider here the case in which function f{x) (0 < x < 1) is 
given by its discrete values, fi = /(A/2 + A • (i — 1)), (i = 1, 2, • • • , 2''), where 
A • 2"" = 1. By the use of column matrices, 



F = 



h 



A,^ 



2 



J =0,1,' 



the wavelet coefHcients can be written as 

oo = 2~iLiL2---LrF, 
A,-i = 2-^HjL.j+i---LrF, j = l,2,---,r-l, 
Ar-i = 2-^HrF. (5) 

where Lj and Hj are 2-'"^ x 2^ matrices. The «— th row of Lj is expressed as, 

( • • • l2[i-l) + l ?2(i-l)+2 • • • '2(i-l)+Ar • • • ) 

= 4= ( • • • Co ci • • • cjv-i • • • ) . 

If 2(i - 1) + A: > 2^ element /2(i-i)+fc is added to the (2(i - 1) + fc - 2^) -th 
element in each row. Matrix Hj is obtained, if Cfc (fc = 0, 1, • • • , A^ — 1) in Lj 
is replaced by {—l)''~^^CN-i-k- 

Matrices Lj and Hj satisfy the following conditions (coming from Eq. (Q)), 

Lj *Lj = Hj *Hj = I, Lj *Hj = Hj *Lj == 0, (6) 

where *ij denotes the transpose of matrix Lj, and / is the 2^^-'^— th order unit 
matrix. 

From Eqs. (|^) and (H), one gets the inverse wavelet transform: 

r-l 

F = F''' + J2p^'\ 

3=0 

^F"^ = 2^ aoLiL2---Lr, 
*F(-''-i) ^ 25 *A,_iiJjL,+i---L,, j = l,2,---,r-l, 

*F('-1) = 25 'Ar-lHr, (7) 



Finally, Eqs. (^ and (|^) lead to the following identity; 



FF2-^ = 2-^/J^a^ + ^i?,-, 

2^-1 

i?, = *A^=^4„ (8) 

fc=0 

0. 



where Ej denotes the j-th level wavelet spectrum 
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3 Daubechies' mother wavelets 

We prepare now Daubechies' mother wavelets, D4,Dq, Ds, and -Dio, cf. Fig. 
1, to be used for wavelet analysis below. If results are independent on these 
mother wavelets, one can say that wavelet spectra bear some physical /or 
stochastic meaning to be further explored. 



Fig. 1. Daubechies' mother wavelets, D4,Dq,Ds, and Dio- 

4 Wavelet spectra of the data 

We analyse now by the wavelet transform LrHuH the pseudo-rapidity distribu- 
tions of the three JACEE events and the TLS event. 

1) In the case of Ca-C and Si-AgBr JACEE events (cf. Fig. 2) both 
original data and subtracted data ( original data - its background (BG)d) are 
used. The corresponding wavelet spectra are shown in Fig. 2. 



Fig. 2. Ca-C and Si-AgBr JACEE events and their wavelet spectra. Results 

obtained by the Poisson random number distributions are shown with l^Ej) + ctj 

and {Ej) — gj. Notice that "(...)" means an assemble average. 



We shall propose now how to distinguish a prominent peak from apparent 
ones. To this end, we use the wavelet coefficients satisfying the following 
condition, 

\ajk- ti\> acr, (9) 

where /i is the average and a is the standard deviation of the wavelet coeffi- 
cients. If a is chosen as 1.96/or 2.58, the absolute values of the coefficients 
become large. Our results show that peak present in Si-AgBr event does not 
depend on the Daubechies' mother wavelets D4,Dq,Ds, and Diq. Therefore 
we can examine further its physical meaning. On the other hand, the peak in 
Ca-C event depends on the Daubechies' mother wavelets. Therefore we can 
say that it is an apparent one. 



Fig. 3. Pseudo-rapidity distribution constructed by larger coefficients. 
2) The TLS event is a photon distribution. Using the same procedure, we 
obtain its wavelet spectra shown in Fig. 4. At present it is difficult to elucidate 
physical meaning of their oscillatory behaviour. 



Fig. 4. Data of TLS event and its wavelet spectra. 



3) In this case of JACEE event G27, Arizona groupEa and Bjorkenllj have 
conjectured that this event probably contains some signal of the disoriented 
chiral condensate (DCC). The ratio r — no/{no + nc/i)tJ is shown in Fig. 5. 
Using the same procedure as mentioned, we obtain the wavelet spectra shown 
in Fig. 5. 



Fig. 5. JACEE G27 event and its wavelet spectra. 
It is found that our result is different from that of Arizona groupEj (Notice 
that our method produces for the random distribution an increasing behaviour 
of wavelet spectrum when j increases.) 

5 Concluding remarks 

We have prepared the Daubechies' mother wavelets, D4, Dq, Dg, and Diq, and 
using them we have analysed high multiplicity events. It is found that our 
wavelet spectra obtained from three JACEE events and TLS event show dif- 
ferent behaviour from that of Poisson random number distributionQ. We con- 
clude therefore that one can use the wavelet analysis to gain some physical 
information/or stochastic properties from the high multiplicity distributions. 
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Fig.2(c) 
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Fig.3(a) 
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Fig.3(b) 
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Fig.3(c) 
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Fig.3(d) 
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